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ABSTRACT 

Context. The optical surface brightness of dark nebulae is mainly due to scattering of integrated starlight by classical dust grains. It 
contains information on the impinging interstellar radiation field, cloud structure, and grain scattering properties. We have obtained 
spectra of the scattered light from 3500 to 9000 A in two globules, the Thumbprint Nebula and DC 303.8-14.2. 
Aims. We use observations of the scattered light to study the impinging integrated starlight spectrum as well as the scattered Ho- and 
other line emissions from all over the sky. We search also for the presence of other than scattered light in the two globules. 
Methods. We obtained long-slit spectra encompassing the whole globule plus adjacent sky in a one-slit setting, thus enabling effi- 
cient elimination of airglow and other foreground sky components. We calculated synthetic integrated starlight spectra for the solar 
neighbourhood using HlPPARCOS-based stellar distributions and the spectral library of Pickles. 

Results. Spectra are presented separately for the bright rims and dark cores of the globules. The continuum spectral energy distri- 
butions and absorption line spectra can be well modelled with the synthetic integrated starlight spectra. Emission lines of Hq'+[Nii], 
Hy6, and [S ii] are detected and are interpreted in terms of scattered light plus an in situ warm ionized medium component behind the 
globules. We detected an excess of emission over the wavelength range 5200-8000 A in DC 303.8-14.2 but the nature of this emission 
remains open. 

Key words. ISM: clouds - ISM: lines and bands - dust, extinction - solar neighborhood - Radiative transfer - Methods: observational 
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1. Introduction 
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.Since the pioneering work by IStruve & Elvevl(ll936h . the sur- 
face brightness of dark nebulae has been recognised as a tool for 
studying the scattering properties of interstellar grains. Broad- 
band photometric studies have been used to constrain the grain 
'albedo and sca ttering phase function in the optical ( Mattilaj 
[1970al : iFitzgerald, Stephens. & Witt 1976). ne ar-IR (Lehtinen 
.& Mattila 1996) and UV (iHaikala et alJll995ir Because of the 
faint surface brightness levels involved, intermediate-resolution 
.photometric or spectrophotom etric studies have been carried out 
■only in a few cases so far (Mattila 1979; Laureijs, Mattila, & 
"Schnur Il987h . In these studies the impinging interstellar ra- 
diation field (ISRF) has been assumed to be known in the so- 
■lar neighbourhood. Broad band observational determinations 
of the optical ISRF include Elsasser & Haug ( 1960) and Tolled 
(119811) . while som e frequently used IS RF models have been pub - 
lished bv lHabind (11968) and Mathis. Mezger. & Panagial (Il983h . 
A higher resolution spectral energy distribution(SED) and line 
spectrum of the optical ISRF has been calculated by 'Mattila' 
( 1980a,b) using population synthesis modelling of the integrated 
starlight. 

Higher spectral resolution observations of the surface bright- 
ness of dark nebulae are desirable for at least three reasons: (i) to 



* Based on observations collected at the European Southern Obser- 
vatory, Chile, under programme ESO No. 073.C-0239(A) 
** Closed down on December 31, 2009 



study the SED of the ISRF in more detail with the aim of justify- 
ing the population synthesis modelling; (ii) to study the wave- 
length dependence of the dust-scattering properties in greater 
detail; and (iii) to find out whether other mechanisms, such as 
photoluminescence, contribute to the surface brightnes in addi- 
tion to scattering. In this paper we are concerned with items (i) 
and (iii). 

The observed optical/near-IR scattered light spectrum of a 
dust cloud, externally illuminated by the Galactic ISRF, can be 
used to empirically determine the spectrum of the ISRF. Indeed, 
the clouds can be used as a kind of ISRF remote sensing measur- 
ing devices or modified Ulbricht spheres (Ulbricht 1920), placed 
in differ ent locations of the solar neighbourhood. Recentlv, Mat- 
tila et al. (I201lh have published the spectrum (3500-6000 A) of 
the high galactic latitude dark cloud LI 642 and found a good 
agreement with the integrated starlight spectrum synthesized for 
the solar neighbourhood. Brandt & Draine (2012) have deter- 
mined the average scattered light spectrum of diffuse dust by 
using 92000 blank-sky-spectra from the Sloan Digital Sky Sur- 
vey. 

Spectroscopic observations are also needed to check whether 
the extended red emission (ERE), detected so far in reflection 
nebulae, planetary nebulae, HII regions, galaxies, and the dif- 
fuse interstellar medium is present also in dark nebulae. A previ- 
ous spectro photo metric measurement of t he dark nebula L 1780 
by Mattila (I1979h has been interpreted bv lChlewicki & LaureiisI 
(il987 ) to show evidence for ERE. Because of the low spectral 
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resolution (R ~ 30) and relatively narrow wavelength coverage 
(5000 to 7500 A) of this measurement, there is obvious need for 
additional, wider band and higher resolution spectral measure- 
ments of the surface brightness of dark nebulae. This need was 
emphasised also in the compilation and comprehensive discus- 
sion of the existing ERE data by Smith & Witt (2002), where 
L1780 had an important place. 

Based on obse rvations made with the ESO Multi-Mode In- 
strument (EMMI) (iDekkeret al.lll986h at the 3.5-m NTT at ESO 
La Silla Observatory, we present spectra from 3500 to 9000 A 
of the scattered light in two globules. Our observations extend 
the number of dark clouds observed for scattered light spectrum 
from one to three and the wavelength coverage up to 9000 A. 

The two globules of this study, the Thumbprint Nebula 
(DC302.6-15.9) and DC303. 8-14.2 (to be called TPN and 
DC 303 for brief), are located in the Chamaeleon III molecu- 
lar cloud area. For the description of the globules and a review 
of previous optical, infrared, and radio molecular line observa- 
tions see |Fitzgerald ( 1974), Fitzgerald, Stephens, & Witt ( 1976), 
iLehtinen et al. (1995), Lehtinen & Mattila (1996), Lehtinen et 
al. ( 1998) and Kainulain en et aL (2007). The TPN has a regular 
elliptical shape with a bright rim - dark core structure. Its appar- 
ent size is ~ 5 X 3 arcmin. DC 303 is somewhat larger, ~ 10 x 5 
arcmin, less regular in shape but with a similarly pronounced 
bright rim - dark core structure as TPN. The density structure 
of T PN and DC 303 has been studied using deep JHK photom- 
etry (iKainulainen et aLll2007h . DC 303 has a dense centrally 
peaked core with extinction Av > 50 mag, while the TPN has a 
more transparent core with a peak extinction of Ay » 13 mag. 
DC 303 has an optically invisible embedded far-infrared source 
near its centre, a Class (3/1 you ng stellar object. The TPN is star- 
less (IKainulainen et al]l2007h . As a distance est imate for bot h 
globules we use the value for Cha III, 217+12 pc (lKnudell2010l) . 

The surface brightnesses in this paper is given in units of 
10"'' ergs cm"^ s ' sr"' A"' , for which we use the abbrevia- 
tion cgs. For the integrated spectral line intensities we use the 



TPN 



abbreviation 1 cgsA^ 



2. Observations 



:10 ^ ergs cm ^ s ' 



sr 



The o bservations were made with the EMMI dDekker et al.l 
[1986') at the New Technology Telescope (NTT) on the La Silla 
observatory during the nights of 17-19 May 2004. The spectro- 
scopic observing modes used were BLMD (blue medium dis- 
persion spectroscopy) in the blue arm of the instrument and 
RILD (red imaging and low dispersion spectroscopy) in the red 
arm. The detector in the red arm consists of two CCD chips 
(2048x4096 pixels) with a gap of 7.5"We used 2x2 binning, 
giving a pixel scale of 0.33"/pixel. The detector in the blue arm 
is a single 1024x1024 pixel chip that was binned by 2 x 2 pix- 
els, giving a scale of 0.72"/pixel. In the BLMD mode we used 
grating #4 as the dispersive element. In the RILD mode we used 
grism #3 as the dispersive element. The slit width was set to 
5". The slit length was ~6' and ~8' in the blue and red arm of 
EMMI, respectively. Details of the instrument setup are listed in 
Table[T] The integration time of a single spectrum is 30 minutes. 
For the intensity calibration of spectra we measured several spec- 
trophotometric standard stars each night through a 5" slit. The 
wavelength calibration of spectra is based on arc lamp frames. 

The slit of EMMI was positioned through the centre of the 
rather spherical clouds in such a way that the slit included the 
dark core, bright rim(s), and pure background sky (see Fig. [TJ. 
The diameter of TPN, ~5', is such that both of the bright rims 
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Fig. 1. Slit positions superimposed on Digitized Sky Survey red im- 
ages of the TPN and DC 303. The length of the slit, 8', is that of the slit 
in the red arm of the spectrometer. The position angle (East of North) 
of the sUt is 0° and 168° for TPN and DC 303, respectively. 



fit inside the slit. Because the diameter of DC 303 is larger than 
the length of the slit, we measured only one bright rim, the dark 
core, and the background sky. Digitized sky survey images were 
used to set the centre position and orientation angle of the slit in 
a way that minimises the number of stars within the slit. 

In this project, we made spectroscopic long-slit observations 
of the weak surface brightness of a dust cloud that extends over 
a considerable length of the slit. Furthermore, we used blank sky 
areas in the science frames to determine the spectrum of the night 
sky emission, which was then subtracted from the object spec- 
trum. For a proper calibration of these data, two items require 
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Table 1. Instrument setup of the EMMI instrument for spectroscopic observations. The spectral resolution, measured from arc-lamp spectra, is 
defined in Sect. 3.3. The wavelength range given is the range used for making the spectrum. 



Mode 


Dispersive 


Dispersion^ 


Spectral 


Wavelength 




element 


(A/pixel) 


resolution (A) 


range (A) 


BLMD" 


grating #4 


3.68 


44 


3381-5119 


RILD'= 


grism #3 


2.86 


25 


4302-9012 



Notes. For 2x2 pixel binning used in the observations Blue medium dispersion spectroscopy Red imaging and low dispersion spectroscopy 



special attention. Firstly, we must be able to make a good flat- 
field correction for the large-scale gradients over the frame. The 
best flat-field correction, both for imaging and spectroscopy, is 
obtained when the illumination of the instrument and the spec- 
trum of the radiation entering the instrument are the same be- 
tween flat-field calibration images and science images. We ful- 
filled these conditions by taking separate flat-field images during 
the nighttime, pointing the telescope to a position near the glob- 
ules that is free of gradients in dust emission, based on 100//m 
IRAS images. Secondly, we must be able to remove the geomet- 
rical distortions in the dispersion axis of the spectra, to ensure 
that the dispersion axis is along a row or column over the whole 
frame. Otherwise, subtracting the night sky emission spectrum 
produces artefacts on the object spectrum, especially at wave- 
lengths where the gradient in the night sky spectrum is at steep- 
est. We performed this geometrical correction by using the arc 
lamp calibration spectra of the EMMI instrument. 

3. Data reduction 

Data reduction wa s done usin g IRAlQ (Image Reduction and 
Analysis Facility) ( lTodvlll993l) . The details of the data reduc- 
tion are the following. 

i) Bias subtraction. We used separate bias images to determine 
the bias level in the science frames. At the blue arm the master 
bias image is flat without structure, so the bias value was taken 
to be a scalar and subtracted. At the red arm the master bias im- 
age shows a gradient across the frame, therefore a 2D bias image 
was subtracted. 

ii) Flat-field correction. To make the flat-field correction for the 
large-scale gradients over the images, we used the nightsky flat- 
fields as discussed in Sect. 2. These flat-field images are weakly 
illuminated and therefore do not have a sufficiently high signal- 
to-noise ratio. Therefore, we first replaced every line of the flat 
field image in the spatial direction with a Chebyshev function fit 
of the line, and discarded any stars in the fit. Each line was then 
scaled to a mean value of unity. We divided the science frames 
with each individual flat-field image and used the flatness of the 
background sky outside the globules to decide which flat field 
image gives the best correction. To flat-field-correct for pixel- 
to-pixel sensitivity variations we used the dome flat images. We 
first took a mean of all dome flat images, then smoothed the 
mean image with a Gaussian function, and then divided the mean 
image with its smoothed version. This provided a flat field image 
that only includes the pixel-to-pixel sensitivity variations. 

We coiTected for geometrical distortion of the spectra in the 
following way in IRAF: i) on a 2D arc-lam spectrum the emis- 
sion features along a single dispersion line are identified with 

' IRAF is distributed by the National Optical Astronomy Observato- 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci- 
ence Foundation 



the identify task of IRAF; ii) the emission features at other 
dispersion lines are re-identified with the reidentify task; iii) 
the wavelengths of the identified features as a function of pixel 
position are fitted with a two-dimensional function using the 
fitcoords task; iv) the geometrical correction is made with the 
transform task, after which the wavelength is a linear func- 
tion along one axis (dispersion is constant). Despite this cor- 
rection, subtracting the background sky emission produces arte- 
facts on spectra at the edges of the brightest airglow emission 
line profiles where the intensity gradient is steepest. Therefore, 
we interpolated the following wavelength ranges of the spectra: 
5544 A-5601 A, 5909 A-5918 A, and 6277 A-6322 A. 



3.1. Atmospheric extinction 

The observed spectra have to be corrected for atmospheric 
absorption lines. To derive a spectrum of telluric absorp- 
tion we used the spectrum of the standard star LTT4364 ob- 
served bv us. The spectral class of LTT4364 is D06.4 (Sion 
et al. 20091) and it is free of significant intrinsic absorption 
features in the wavelength range 6000-9000 A. At the spec- 
tral resolution of our observations, about 44 A, the individ- 
ual, narrow absorption lines blend into wide absorption bands 
that have wide and weak wings (see Fig. |2]i. Thus it is dif- 
ficult to determine where the continuum (areas not affected 
by telluric absorption) in the spectrum of LTT4364 is lo- 
cated. To determine the wavelength ranges of the continuum 
we used the tabulated atmospheric transmission spectrum pub- 
lished by NOAO (National Optical Astronomy Observatories) 
(ftp://ftp.noao.edu/catalogs/atmospheric_transmission). Firstly, 
we convolved the NOAO spectrum to the resolution of our ob- 
servations. A comparison between the convolved NOAO spec- 
trum and the LTT4364 spectrum verifies that the absorption fea- 
tures seen in the LTT4364 spectrum are all caused by telluric 
absorption. Secondly, we defined the continuum as wavelength 
ranges where atmospheric transmission in the convolved NOAO 
spectrum is higher than 99%. Then, we fitted the observed spec- 
trum of LTT4362 with a Chebyshev polynomial, excluding the 
data points that are outside the continuum, as defined above. 
The observed spectrum of LTT4364 and the continuum fit are 
shown in Fig. |2^. By dividing the spectrum of LTT4364 with 
the fit, we obtain the normalised telluric transmission. As the fi- 
nal steps, we i) applied a Savitzky-Golay smoothing filter to the 
derived transmission spectrum (reducing noise while retaining 
the dynamic range of the variations); ii) defined the transmis- 
sion to have a value of unity in the continuum regions as defined 
above; and iii) scaled the transmission to transmission at air- 
mass of unity by raising the transmission to a power of (l/X) 
where X is the airmass of LTT4364 during our observation. The 
final transmission spectrum is shown in Fig.|2}3. At wavelengths 
longer than about 8900 A the derived transmission is extrapola- 



Article number, page 3 ofl20l 




0.5 h I i 

oa I ............................. .1 

6000 6500 7000 7500 8000 8500 9000 
Wavelength [A] 

Fig. 2. Derivation of atmospheric transmission. The observed spec- 
trum of the standard star LTT4362 (histogram), overplotted by a Cheby- 
shev fit (solid line) to the continuum of the spectrum. The regions that 
are defined to represent the continuum are shown as intermittent lines 
above the fit (a). The derived atmospheric transmission spectrum at air- 
mass of unity (b). The molecular bands of the spectrum are identified. 



tion because of the absence of a continuum region at the long 
wavelength side of the absorption feature. 

Each individual globule spectrum was corrected for telluric 
absorption by dividing by the transmission spectrum (Fig. |2j5) 
raised to the power of X, X being the airmass of the globule. 
Correcting the deep absorption feature at ~ 7600 A leaves no- 
ticeable residuals in the spectra, therefore we interpolated the 
wavelength range 7568 A-7652 A after the telluric absorption 
correction. 

3.2. Extracting the spectra 

The main products to be derived from our spectroscopic data are 
the spectra of scattered light at the position of the bright rims and 
the dark cores. Because the optical depth through a dust cloud 
depends on the wavelength, the shape of the surface brightness 
profile changes accordingly (see Fig.|4|i. Thus, one should deter- 
mine the extent of the bright rim and the dark core in the spatial 
direction in a way that is dependent of the wavelength. 

In the red arm we used the following method. Firstly, to 
obtain a sufficient signal-to-noise ratio, we rebinned the surface 
brightness profiles over wavelength ranges of ~33 A (1 1 pixels). 
The rebinned profiles were fitted with a Chebyshev function in 
spatial direction, and the pixel locations of the maxima of the two 
bright rims and the minimum of the dark core were noted. The 
bright rim was defined as a window over those pixels in spatial 
direction where the surface brightness exceeds the maximum of 
the bright rim by 0.8 times, and the dark core was defined as a 
window where the surface brightness is less than 1.2 times the 
minimum of the dark core. The derived extents of the rims and 
the core as a function of wavelength were fitted with a low-order 
polynomial, and the fit was used to interpolate the windows at 



every wavelength of the full-resolution spectra. The windows of 
the bright rims and dark core defined in this way are shown in 
Fig.|4]for a single spectrum of TPN. 

In the blue arm the signal-to-noise ratio in individual spectra 
is so low that this method is not applicable. Therefore, we used 
constant ranges in spatial direction to specify the extent of the 
bright rim and the dark core. 

We then took the mean of the pixel values at each wavelength 
over the spatial windows as defined above. The blank sky at the 
northern side of the globules was used to determine the inten- 
sity level of the background sky at each wavelength, which was 
then subtracted. We defined the intensity AI as an intensity of 
diffuse light of the globule relative to the intensity of the back- 
ground sky. To improve the signal-to-noise ratio, we rebinned 
the spectra into 20 A wide bins. 

3.3. Wavelength resolution 

The point spread function (PSF) of the spectra in dispersion di- 
rection resembles a square-box profile, as determined from emis- 
sion line profiles in arc-lamp measurements. Therefore we give 
the wavelength resolution as effective width, determined as the 
area under an emission line profile divided by the maximum of 
the profile. In the red arm, the resolution changes from about 
42 A to about 45 A from the short- to the long-wavelength edge 
of the spectra. The mean value of the resolution over 13 separate 
emission lines of an arc-lamp spectrum is 44± 1 A. In the blue 
arm, the resolution over four separate emission lines is 25±1 A, 
with no dependence on wavelength. 

3.4. Error analysis 

Most of the total signal measured by us is nightsky emis- 
sion. Consequently, the total number of detected photons varies 
strongly as a function of wavelength (see Fig.|4]i. The maximum 
surface brightness of the diffuse light from the globules is *:1/10 
of the brightness of the pseudo-continuum of nightsky emission 
(that is, outside the emission peaks). In the red arm of EMMI, 
the total number of detected electrons in a single exposure varies 
between ~300-2500 e" per pixel, depending on the strength of 
atmospheric emission. For the pseudo-continuum the intensity 
is typically about 400 e" per pixel. The signal-to-noise ratio for 
diffuse light at the bright rims is thus about two per pixel on a 
single exposure. When combining all individual exposures, we 
have typically ~2000 pixels on spatial dkection at each wave- 
length over which to calculate the mean value of the diffuse light 
intensity at the bright rims. The standard deviation of the dif- 
fuse light spectrum of the bright rim caused by photon noise is 
thus ~ 1 % of the intensity of scattered light. This deviation can 
be considered negligible compared to other noise sources. The 
other sources of eiTor that we considered are related to photomet- 
ric calibration and imperfect flat-field correction and background 
subtraction, producing an estimated total uncertainty of 25%. 

4. Results 

4. 1 . Surface brightness profiles 

Figure |3] shows examples of the surface brightness profiles for 
TPN and DC 303 as measured along the slit, after subtracting 
the background sky emission. The region at the leftmost part 
of the plots represents background sky which has been used to 
determine the zero level of intensity at each wavelength. The 
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Table 2. Detected and searched lines and bands in the spectra of scattered light. Results are given separately for rim and core spectrum. Columns 
3-6: line is seen in emission ('em') or absorption ('abs') relative to the background sky, or not detected at all (blank). Columns 7 and 8: strong 
atmospheric emission or telluric absorption make the status of the line uncertain. 



Line or Wavelength TPN rim TPN core DC 303 rim DC 303 core Atmosph. emission Atmosph. absorption 

band [A] contamination contamination 



1 


2 


3 


4 


5 


6 


7 


G-band+H-y 


4300+4341 


abs 


abs 


abs 






H/3 


4861 


abs 


abs 


abs 






Mgl + MgH 


5180 


abs 


abs 


abs 






NaDl 


5896 


abs 




abs 


abs 


V 


NaD2 


5890 


abs 




abs 


abs 


V 


TiO 


[7050-7450] 


abs 




abs 




V 


Ha+[Nii] 


6563+6548/6583 


em 


abs 


em 


abs 




[Sii] 


6716/6731 




abs 


em 


abs 




TiO 


[6150-6300] 


abs 


abs 


abs 


abs 


V 


Can 


8498+8542 


abs 


abs 


abs 


abs 





Can 8662 V 



profiles are for a band containing the Ha and [N n] lines, together 
with the underlying continuum, and for an adjacent continuum 
band (mean value over the wavelength ranges 6473-6534 A and 
6615-6676 A). Also shown is the ratio of these two bands. The 
surface brightness profiles were scaled to the maximum intensity 
of the band containing the Ha and [N n] lines. 

4.2. Spectra of the bright rims and dark cores 

The spectra of scattered light are shown in Figs.|5]and|6]for TPN 
and DC 303. Panels a show the spectra for rim and core, while 
panels b show the ratio of rim to core. In the following we give 
a description of the observed properties of the lines and bands 
that were detected or searched for in our spectra. 
Ha+[Nii]: The most prominent line is Ha+[Nii]. It is seen in 
emission in the rim of both clouds. In the core spectra, the Ha 
line is seen in absorption for both clouds. With our spectral 
resolution the Ha line includes the [Nii] lines at 6548 A and 
6583 A, which were detected in emission in the diffuse galactic 
light by Brandt & Draine (20i3)- 

[Sii]: The second-most prominent line is [Sii], including two 
lines at 6716 A and 6731 A, which are unresolved with our 
resolution. In TPN the [S ii] line is absent in the rim spectrum 
and negative in the core spectrum. In DC 303 the [S n] line is 
positive in the rim spectrum and negative in the core spectrum. 
Hy + G band: These two features are blended in our spectral 
resolution. In TPN these lines are not prominent, but the ratio 
spectrum shows a negative peak close to the G band. In DC 303 
the lines can be seen in absorption in the rim spectrum, but are 
absent from the core spectrum. The ratio spectrum shows a 
wide negative peak ranging from ~4250^500 A that cannot be 
explained solely with Hy and the G band. The Fraunhofer G 
band is due to absorption by the CH molecule, with contribution 
from other molecular bands and metals. 

H/?: This line is seen in absorption in the rim and core spectra of 
TPN and in absorption in the rim spectrum of DC 303, while it 
is absent from the core spectrum of DC 303. 
Mgl + MgH: This band is seen in absorption in all spectra 
except the core spectrum of DC 303, from which it is absent. 
Na D1+D2: The unresolved combination of these lines is in 
absorption in all spectra except the core spectrum of the TPN. 
TiO bands: These bands are in absorption in all spectra. 
Can triplet: The blend of the 8498 and 8542 A fines of this 



triplet is in absorption in all spectra. The third fine at 8662 A is 
disturbed by a strong airglow band. 

A compilation of the spectral features is given in Table|2] In 
addition, the spectra include features that are artifical because 
they are only one channel wide. The emission-like feature at 
~6900 A is caused by the edge of the telluric O2 absorption band. 

4.3. Ratio ofbrigth rim and dark core spectra 

Because the main noise contributor is the same in the spectra of 
the bright rim and the core, the noise between them is correlated. 
Therefore, the noise is expected to be reduced for the ratio of the 
bright rim and the core spectra. Furthermore, for the TPN the 
rims are located symmetrically around the core, thus any imper- 
fections in the correction of geometrical distortions (see Sect. 3) 
are largely cancelled out in the ratio. Because the intensity of 
the core spectrum crosses zero, we derived the ratio in the fol- 
lowing way: i) We fitted the spectra with a Chebyshev function, 
separately for wavelengths below and above the 4000 A break, 
ii) We subtracted the fitted function from the spectra and scaled 
the spectra so that the mean value of the rim and core spectra is 
unity, iii) Finally, we divided the rim and core spectra. 

5. Modelling scattered light from globule and 
adjacent sky 

5.1. Qualitative aspects 

To explain qualitatively the behaviour of scattered light at the 
bright rim and core of the TPN and DC 303 the following facts 
are to be realised. 

i) The sky adjacent to the globules has a substantial dust col- 
umn density. Accordingly, the background sky has continuum 
and line emission caused by scattering of the interstellar radia- 
tion field, plus a signal from the in situ gas emission. The extinc- 
tion for the skv areas was determined using 2MASS (Skrutskie 
et al. 2006) JHK photo metrv and the N ICER method of Lom- 
bardi & Alves (2Q0D: Ay(sky) ^ 1 mag for both the TPN and 
DC 303. The extinction and scattered light in front and behind 
of the globules have to be considered separately. The foreground 
extinction layer can be assumed to be homogeneous over the slit 
length. Then, its contribution to the scattered light difference is 
cancelled out but influences the measured signal via extinction. 
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Fig. 3. Surface brightness profiles across the globules in a band con- 
taining the Ha and [N n] lines, together with the underlying continuum 
(solid line) and the adjacent continuum (dashed line) as well as their 
ratio (histogram). Panel a is for TPN, panel b for DC 303. The narrow 
peaks in the profiles are stars. The ratio is shown only for the region 
where the continuum intensity is not close to zero. The outer, northern 
boundary of DC 303 is on the left-hand side of the bright rim at ~0.7'. 



The location of the globules with respect to the extinction layer 
is not known. Because most of the dust in this direction is lo- 
cated at the distance of the Chamaeleon cloud complex dKnud^ 
I2OIO) . we assumed that the globules are located in the middle 
of this layer, with equal fore- and background extinctions of 
Ay = 0.5 mag. 

ii) At the position of the bright rim the optical depth through 
the globule, t(/1), has a value of ~2 mag. The conditions are 
therefore ideal for the scattering process. At the very shortest 
wavelengths of our observations the intensity of the rims above 
sky background approaches zero. This is mainly because the 
scattered light intensity at the sky position is increasing towards 
shorter wavelenghts proportionally to 1 - e^'^^'^\ while the brigth 
rim signal is saturated and remains proportional to the ISL spec- 
trum. A decreasing albedo of dust particles may also contribute. 

iii) At the position of the cores the optical depth through the 
TPN and DC 303 is much higher than unity, t{A) > 10, at all 
wavelengths covered by our spectra. The scattering process is 
therefore influenced by multiple scattering (and absorption) and 
the intensity of the scattered light is consequently reduced. Fur- 
thermore, the light from behind the core is blocked. Thus, the 
surface brightness of the core relative to the sky becomes nega- 



tive below a threshold wavelength. For the TPN this happens at 
A ~ 4300 A, for DC 303 at A ~ 5900 A. 

iv) The scattering geometry is different for the bright rims 
and the dark core. Interstellar grains are known to have a 
strongly forward-directed scattering function in the optical band, 
with the commo nly used Henvev-Greenstein parameter (Henyey 
& Gre enstein 1941) in the range g ^ 0.7 - 0.8 (Mattila' 1970b'; 
IWitt. OUveri . & Schild 1990: Drain3l2003h . This means that for 
the sky areas as well as for the bright rims the illumination is 
coming predominantly from behind the globule. On the other 
hand , the increasing optical thicknes has been shown by Mattil^ 
(1970b) to modify the "scattering function of the cloud" from 
forward directed towards a more isotropic one. Thus, for the 
dark cores that are optically thick (Ay > 10™) the illumination is 
coming more uniformly from all over the sky or even from the 
backward direction. 



5.2. Nature of the Ha emission from the globules: scattered 
light 

If the Ha and other line emissions from the globule rims were 
caused by an ionised layer on the surface of the globule, one 
would expect to see a narrow intensity peak at the outer surface. 
On the other hand, if the emission is mainly scattered light, its 
distribution should be similar to the distribution of the adjacent 
continuum. We show in Fig. [3] the surface brightness profiles 
of a band containing the Ha and [N 11] Unes, together with the 
underlying continuum, and a band of the adjacent continuum. 
Also shown is the ratio of these two bands accross the bright 
rims in each of the globules. Clearly, there is no hint of a thin 
ionised gas layer at the outer boundaries of the globules. The 
profiles over the bright rims are consistent with the idea that the 
Ha and [N 11] line emission in the bright rims is also scattered 
light. 

5.3. Modelling the continuum spectrum with ISL 

Because the observed spectra are differences between the glob- 
ule and the nearby sky, the large foreground components of the 
sky brightness, the zodiacal light, airglow, and tropospheric scat- 
tered light are eliminated. 

As explained in the appendix, we calculated for the solar 
neighbourhood synthetic integrated starlight (ISL) spectra at dif- 
ferent galactic latitudes and for locations at the Galactic plane 
(z = 0) and at the height z = 50 pc corresponding to the location 
of the two globules. 

The scattered light spectra of the globule rim and core are 
modelled as 



rUA) = IisiWCssooe- 



■' (e/)[fl(i)-fl(5500A)] -Tv(fg)aW 



(1) 



whe re IrsiiX) is the ISL intensity, a(A) is the interstellar extinc- 
tion (CarSCci^tonZ&M^hi3i98^ = 3.1 normahsed 
to 1 at 5500 A, Tv(fg) the foreground extinction and T^y^{cl) the 
effective extinction of the cloud (rim or core) at 5500 A in op- 
tical depth units. C5500 is a scaling factor to adjust the intensity 
level to the observed spectrum. 

The scattered light intensity of the background dust layer at 
the sky position is modelled as 



jsky.bg. 



(A) = /;L^'''^(5500A) 



//sz.(5500A)(l -e-^''(''9))' 



-Tv{fg)a(A) 



(2) 
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where Tv(bg) is the extinction at 5500 A of the background 
dust layer in optical depth units. The observed continuum sur- 
face brightness globule-minus-sky is then given by 
AI.M = 4,{A) - Ci'^'iA). 

The sky background intensity at 5500 A, if^'f'^ (5500 A), has to 
be estimated independently. We used two different methods: 

(i) For an opaque core the scattered light in the ultravio- 
let, A :$ 4000 A, approaches zero. This is expected to be 
the case especially for DC 303, which has a core extinction of 
Av > 50 mag. The observed Al^ca at 4000 A is ~ -8 cgs. 
After correcting for the foreground extinction of 0.7 mag at 

4000 A we obtain /J*f'''*(4000 A)=15 cgs, and the same value 
for 4^^^(5500 A). 

(ii) We also used our surface photometry results for the high 
galactic latitude dark cloud L 1642 (/ - 210 deg, b ^ -37 d eg) 
(iLaureiis. Mattila. & Schnujll987l : lMattila et al.ll201 lLl20T2h . In 
contrast to the TPN and DC 303, the surrounding sky of L 1642 
is almost free of dust. Ay ~ 0.15 mag, rlca(5500 A) x 3.4 cgs. 
Plots of IscaiA) vs. A(A), corresponding to the I sea vs. rioo^m 
plots in lLaureijs. Mattila. & Schnud (Il987h . Fig. 7, show a lin- 
ear increase of IscaiA) for A(A) < 1, a flat maximum reached 
at A{A) Si 1.5-2 mag, and a decreasing surface brightness at 
A{A) S 2 mag. We used these curves to derive for L 1642 the 
differences 

A/5ca(L 1642) = /jcfl (maximum) - IscaiAy = 0.5 mag), 
which correspond to the differential measurement "bright rim 
minus sky" for the TPN and DC 303. Comparing A/,ca(L 1642) 
values at the four wavelengths 3850, 4150, 4670, and 5470 A 
with the corresponding bright rim values for the TPN, we find 
an almost constant scaling factor of 1.5 ± 0.1 which, when cor- 
rected for the foreground extinction Ay = 0.5 mag, corresponds 
to 1.9 + 0.13. Using this scaling factor and the value of 7 cgs 
for /s£.fl(5500 A) at Ay - 0.5 mag in L 1642, we obtain for the 
TPN background sky brightness the estimate /^'^a'''^(5500A) - 
1.9 X 7 cgs=13 cgs. 

We show in Figs.|7]and[8]the fitting of the observed spectra 
of the TPN and DC 303 with ISL spectra. For each figure the 
upper panel shows the fit for the bright rim and the lower panel 
for the dark core. The sky background spectrum was determined 
according to the recipe as described above and is fixed to the 
value rf^'f'^ = 15 cgs at 5500 A for both globules. The ISL 
spectrum used for the bright rims is the mean of the ISL spectra 
at 1^1=1 1.5 and 23.6 deg. For the dark cores the ISL spectrum is 
the mean over the sky. The values of C5500 and Ty^^(cl) in Eq. (1) 
are adjusted to optimise the fits to the observed spectra and are 
given in the figure legends. Residuals for the observed-minus- 
model spectrum are shown in the lower parts of the figures. For 
the DC 303 rim and more weakly for the TPN rim an excess 
emission above the scattered light model spectrum is found be- 
tween ~5000 and 8000 A. Its properties and interpretation are 
discussed in Sect. 6.4. 



5.4. Scattered and in situ emission iines 

Scattered emission line radiation from all over the sky con- 
tributes to the rim and core spectra of the globules as well as 
to the spectrum of the sky position in a similar way as the ISL 
continuum light. In addition, there is in situ line emission from 
ionised gas in front and behind the globules. The emission 
lines to be considered here are Ha-i-[Nn]/i6548,6583, HyS, and 
[Sii].16716,6731. 




4000 



6000 
Wavelength [A] 




6000 
Wavelength [A] 



Fig. 7. Fitting of the TPN spectra with ISL model spectrum. The 
upper panel is for the bright rim. The observed spectrum is shown 
as a histogram (thick black line). The upper and the lower continu- 
ous thin line in black show ISL model spectra for the bright rim and 
the sky positions, respectively. Their difference, the red continuous 
line, is the fit. The parameters used for the model spectrum were 
A'/^(cloud) = 1.086//^(cloud) = 0.33 mag and C5ooo=0.66. The 
residuals observation-minus-fit are shown as histogram (in red), shifted 
down by 20 units. Tlie lower panel shows the core. Same as for 
the rim but the parameters used for the model fit were A^^-' (cloud) = 
1.086t/-' (cloud) = 1.4 mag and C5ooo=0.33. 



These lines, after subtracting the integrated starlight model 
spectra, are found to be positive (in emission) or zero in the rim 
and negative (apparently in absorption) in the core spectra, be- 
cause emission lines from the sky position dominate in the latter 
case. We denote with /'^^ the scattered line intensity integrated 
over the emission line profile and with /|„j„„ the intensity of the 
in situ line emission from behind the globule. The in situ as well 
as scattered line emission in front of the globule is assumed to be 
constant over the 6 or 8 arcmin slit length and thus to cancel out 
in the differential measurements. It is meaningful to combine the 
scattered line component of the background sky behind the glob- 
ule, /'^.g(sky), with the rim or core scattered line component be- 
cause it derives its illumination from the same source, the all-sky 
line radiation dominated by the bright HII regions in the Galactic 
plane. On the other hand, the in situ background line emission is 
produced in the local warm ionised medium (WIM), which has 
different physical condit ions. The two components have diff er- 
ent emission line ratios (iMadsen. Revnolds. & Haffnerll2006l) . 
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Fig. 8. Fitting of the DC 303 spectra with ISL model spectrum. 
Same as for the TPN in Fig.|7]but the fit parameters were A^^-' (cloud) = 
1.05 mag and C5ooo=0.60 for the rim (upper panel) 
2.7 mag and C5ooo=0.10 for the 



1.086//^ (cloud) 
and Af^ (cloud) = 
core (lower panel). 



1.086T/f^ (cloud) 



We assumed that the optical depth through the bright rim at 
the wavelength of each line is r(rim)=2 while the core is as- 
sumed to be opaque, T(core)> 10, transmitting practically no 
background light. The optical depth of the foreground dust layer 
is denoted with T(fg) and corresponds to the dust extinction of 
Av-0.5 mag as adopted in Sect. 5.1. 

The scattered component of the line intensity globule-minus- 
sky can be expressed as 



U.scJrim) + I,^Jsky)e 
i^'scaicore) - l[jsky))e 



-T(rim) 



-T(fg) 



-T(fg) 



We furthermore assumed that the equivalent width, EW, for a 
given line, e.g. Ha, is the same for each scattered line compo- 
nent, i.e. rim, core and sky background and, that it is also the 
same for the TPN and DC 303, 



I'icore) 



sea 
jl 



EW X I'Jrim) 
EW X I'.^aicore) 
EWxr^Jsky), 



where the continuum intensity Z^^^, is at the wavelength of the 
emission line. Thus we can write 



EW X Micairim) 



Al[^Jcore) - EW X A/j^/core). 

The observed line intensity (scattered -i- in situ) difference is 
given by 



A/'(r/m) = £W X A/L(r/m) + 4„,„(e 
M'icore) ^ EW X MlAcore) - I'e 



-T(nm) _ 1 N 
' „-T(/9) 



-r(Ig) 



In these equations the diffential continuum intensities (rim mi- 
nus sky and core minus sky), AF^^^{rim) and AF^^^icore), rep- 
resent the scattered light values fitted to the observations (red 
curves in Figs. Q and [8]) and they were previously corrected for 
the foreground extinction T(fg). 

Given the observed line intensities A/' (rim) and A/'(core) for 
the TPN and DC 303, we solved the group of four equations 
for the three unknowns, EW, /L»«(TPN), and l'i„^,iDC 303). 
The results are given in columns 6 and 7 of Table |3] The three 
scattered line components l[^^{rim)e^'^^J^^\ /'^^(core)e"'"^-'*', and 
l[^Jsky)e^'^^f'-'^ were calculated using the EW and the corre- 
sponding scattered continuum components. They are listed in 
columns 8-10 of Table [3] 



6. Discussion 

6.1. Milky Way ISL spectrum in scattered iigiit 

The SED of the bright rim of a globule is expected to be only 
weakly infuenced by interstellar reddening in the cloud, i.e. that 
it is almost a direct copy of the SED of the impinging ISRF in 
spite of the substantial optical thickness. This is true as long as 
the albedo and scattering function asymmetry parameter g are 
constant over the wavelength range in question. This appears 
to be a good approximation over the range considered here, A - 
3500-9000A dMattilall970bl:lLaureiis. Mattila. & Schnuiil987l: 
lDraindl2003l) . Thus, we expect that the bright rim spectrum is a 
good tracer of the ISL spectrum not only for the Fraunhofer lines 
and spectral discontinuities, but also regarding the shape of the 
SED. 

The SED of the TPN rim comes close to this expecta- 
tion: the best-fit ISL model spectrum (with z = 50 pc) 
is obtained with a reddening correction corresponding to 



A'/-^(cl)=1.0864''(cl)=0.33 mag oxE{B-V) = 0.11 mag. This 
reddening can easily be explained with the slightly decreasing 
albedo, from (3.52 at 5500 A to 0.45 at 3450 A ( Laureiis. Mat- 
tila, & Schnurll987t). 

The SED for DC 303 rim is clearly redder than for the TPN, 
corresponding to A^^(cl)=1.05 mag or E{B - V) - 0.34 mag. 
Since the TPN and DC 303 are exposed to about the same ra- 
diation field, the different rim reddening is probably caused by 
different dust properties. A difference in dust composition is 
also suggested by the presence of a stronger ERE component in 
DC 303 as compared to the marginal detection in the TPN. An- 
other reason could be the different d ensity distributions of the 
two globules (iKainulainen et al.ll2007h . DC 303 being more cen- 
trally peaked than the TPN. A detailed radiative transfer mod- 
elling would be required to study this effect, but this is beyond 
the scope of this paper 

The scattered light intensities for the TPN and DC 303 rims 
are very similar, as can be seen from Figs.|2]and[8]and the scaling 
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factors, C55oo=0.66 and 0.60 for the TPN and DC 303, respec- 
tively. 

The bright rim spectra show the stronger Fraunhofer Hnes 
with very similar strengths as observed in the ISL model spec- 
trum. The 4000 A discontinuity is also pronounced. We regard 
that this, as well as the agreement between the observed and 
modelled SEDs, supports the choice of parameters used in our 
population synthesis model for the ISL. One possible deviation 
may be indicated for the TiO band at ~7100 A. The observed 
band strength appears to be somewhat stronger than predicted, 
but the presence of ERE covering this wavelength range may be 
an explanation as well. This band is sensitive to metal abun- 
dance; thus, adding more metal-rich late-type giants might im- 
prove the fit. 

For the dark cores of the TPN and DC 303, where Ay > 10"', 
the scattered light SED is strongly influenced by multiple scat- 
tering and absorption. Between the subsequent scattering events 
the light experiences wavelength-dependent extinction and the 
scattered light becomes strongly reddened. The reddening law 
as determined for point sources is not valid for diffuse radi- 
ation penetrating through a scattering medium (Mattila 1976). 
For our modelling purpose, however, we formally used the point 
source reddening curve to fit the ISL model SED to the observed 
spectra. As shown in Figs. |7]and|8] good fits can be achieved 
this way. The formal effective extinction values adopted were 
Ay^-'^(cl)=1.086T^''^(cl)=1.4 mag and 2.7 mag, and the scaling 
factors C55oo=0.33 and 0.10 for the TPN and DC 303, respec- 
tively. 

Because the surface brightness at the sky position becomes 
a substantial fraction of, or even exeeds, the core surface bright- 
ness, the scattered Fraunhofer lines become weak in the core- 
minus-sky spectra. On the other hand, because the cores are 
opaque to the background light, their spectra show the emission 
Unes of the in situ gas at sky positions as apparent absorption 
lines. 

6.2. Emission lines: comparison witii other studies 

Equivalent width s of th e scattered emission lines according to 
iBrandt & Draind (1201 2|) are presented in the lower part of Ta- 
ble|3] The agreement with our value for HaH-[Nn]/l6583,/l6548 
is surprisingly good considering the different approaches of the 
two studi es. The fSiil/l6716./l6731 equivalent width of Brandt 
& Draine (I2012h . on the other hand, is significantly larger than 
our value. Two reasons can play a role: (1) The illumination for 
the scattered line radiation is coming from different areas of sky, 
predominantly from behind t he globules in our case a nd more 
uniformly all over the sky i nlBrandt & Draind (l2012h : (2) As 
stated by IBrandt & Draind (1201 2|) there mav be a contribution 
by in situ gas emission correlated with the 100 yum fluctuations 
in their values, which increases the apparent scattered light sig- 
nal. Supporting this, a correlation of Ha with FIR dust emission 
has been claimed bv iLaeache et al.l (0.999) to be evidence of a 
first detection of dust emission by WIM.0 

We give in Table [3] also the Hg+rNiiU 6583,/16548 equiva- 
lent width as measured bv iWitt et al.l (1201 Ol) in the high-latitude 
dark nebula L 1780. Their 60 A wide filter centred at the Ha line 
encompassed also the two [Nii] lines. The [S n]/i6716,/l6731 
equivalent width was derived using line ratio measurements dis- 
cussed in Sect. 6.3 below. The HQ'H-[Nn]/l6583,.16548 equiva- 
lent width is somewhat smaller than found by us from the TPN 



- However, these authors omitted the scattered Ha radiation, which 
may have been the reason for the detected correlation. 



and DC 303 and by IBrandt & Draind (1201 2h from scattering 
by diffuse dust. Part of this difference may be caused by the 
difference in calibration of our spectroscopic and their narrow 
band photometric data. Another reason could be the background 
in situ emission, which was not explicitly handled in the anal- 
ysis of IWitt et all (l2010l) . The equivalent width of 4.5 A for 
[S Ii]/I6716,/I6731 is compatible with our va lue of 3.7 + 1.1 and 
significantly lower than the value of Brandt & Draind (|2012|) . 

The equivalent width determined for the scattered 
HaH-[Nii]/l6583,/l6548 emission in the two globules is 
20.2 A. The mean continuum intensity of the Milky Way 
from Pioneer 10 red photornetry ( with niy < 6.5 mag stars 
included) is 167 SIO© (lTolleilll98lh . corresponding to 173 cgs 
at 6560 A. Thus the HQ-l-[Nn] EW of 20.2 A corresponds to 
3490 cgsA or 14.5 R. Assuming that 26% of this signal is due 
to fNi iU6583,.16548 as found bv .Madsen. Revnolds. & Haffned 
(l2006h for HII regions, we obtain for the mean Ha intensity the 
value of 11.5 R, which agrees reasonably well with the mean 
intens ity of 8.0 R as found from the Ha map of Finkbeine^ 
(l200l . 

In the SHASSA survey (IPaus tad et ani200lh the Ha sky in- 
tensities next to the two globules are 6.4 R (1540 cgs A) for the 
TPN and 7.2 R (1730 cgsA) for DC 303. These values refer to 
the whole line of sight. Our estimates, referring to Ha in situ 
emission from behind the globules, are 56% and 32% of these 
values. This is in accordance with our assumption that the glob- 
ules are embedded half way within the dust and gas layer towards 
the Chamaeleon complex. 



6.3. Emission line ratios 

We list in Table|4]the observed line ratios for the scattered and in 
situ line components. The [Nii]/Hy6 ratio for the in situ case was 
calculated from (Ha+[Nii])/H;6 by assuming that Ha/Hj6=2.87 
corresponding to case B for r= 10000 K and no extinction ( Os- 
terbrock & Ferland 2006j). Likewise, [Nn]/Ha was calculated 
from [N ii]/Hy6 using the same Ha/Hj8 ratio. 

For com parison we give i n column 6 the line ratios ac- 
cording to ferandt & Draine (2012) for the sc a ttered light 
from diffuse dust. Madsen. Revnolds. & Haffneii (l2006h have 
made extensive observations of the line ratios [Nii]/16583/Ha 
and [S n]]/16716/Ha, both for bright HII regions and for high- 
latitude WIM gas. To compare their results with ours, we in- 
cluded both line components of [N ii] and [S n] using the ratios 
[Nii]^6583M6548= 2.97 and [S ii]/l67 1 6//1673 1 = 1 .45 . 

Witt et al. (1201 oh have reported measurements with the Wis- 
consin H-alpha Mapper (WHAM) of the high galactic latitude 
dark nebula L 1780. The line ratios of the excess surface 
brightness averaged over the one-degree WHAM beam were 
[S ii].16731/Ha=0.19 and [S ii]/l6731/Ha=0.16. Using again the 
ratios [Nii]^6583//16548=2.97 and [S n]^67 16/^673 1 = 1. 45 we 
obtain for [S ii]^6716,/16731/(Ha,[Nn]/}6583-i-^6548) the value 
of 0.31 (column 9 of Table |4]l. This is the same value as found 
by us for scattered light in the TPN and DC 303. 

The line ratios [N ii]/Ha and [S n]/Ha differ systematically 
for bright HII r egions and the diffuse high-latitude WIM gas 
(iMadsen. Revnolds. & Haff'neri 12006). Our observed ratio 
[S ii]/l6716,/l6731/(Ha-i-[Nn]/l6583,^6548), is 0.31 for the scat- 
tered and 0. 51 for the in situ component. These values fall in 
between the lMadsen. Revnolds. & Haffne^ (|2006') values for HII 
regions (0.14) and WIM gas (~ 0.7, with substantial variations). 
While the scattered line radiation is dominated by illumination 
from HII regions, a substantial part of it is likely coming from 
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the faint but widely distributed WIM, which causes the enhanced 
line ratio above the HII regions value of 0.14. The value of 
[Brandt & Draine (2012), supposed to be dominated by scattered 
radiation from HII regions, is substantially higher and almost 
the same as our value for the in situ emission. As stated above, 
this may indicate that p art of the [S n] line emission detected by 
lBrandt&Drain3(l2012h is not scattered, but in situ emission. 

Our (Ha+[Nn])/Hj6 ratio is higher for the scattered than for 
the in situ line emission, likely because of the interstellar redden- 
ing of the light from HII regions. The lBrandt & Draind (1201 2f) 
value is consistent with our scattered light ratio. Our [N n]fH/3 
and [Nii]/Hq' ratios have substa ntial scatter. They are consis- 
tent with the HII regions value oflMadsen . Reynolds. & HafFneii 
(2006) but not with the WIM. The Brandt & Draind (1201 2|) value, 
on the other hand, is consistent with the WIM but not with the 
HII regions origin. 

We conclude that the emission line spectra of the bright rims 
and dark cores of the TPN and DC 303 can be modelled as a 
combination of external illumination, leading to scattered light 
from the globules, plus in situ emi ssion by gas behind the glob- 
ules. In th e case of the Ha line, Mattila. Juyela. & Lehtinen 
(l2007h and lLehtinen. Juvela. & Mattilal (l2010l) have previously 
interpreted several observed dark cloud maps in terms of com- 
bined scattered plus in situ emission. Witt et al. (2010) have used 
the low [S iil/Ha ratio observed in L 1780 as supporting evidence 
for the scattered light origin of most of the Ha luminosity of this 
cloud. 



6.4. The excess emission at 5200-8000 A 

The residuals of the observed globule spectra fitted with the syn- 
thetic ISL spectrum are shown in Fig. |9] For the DC 303 rim 
there is an excess emission of ~ 7 cgs over the wavelength range 
5200-8000 A while the TPN rim shows weaker excess emis- 
sion of ~ 3 cgs over 6200-8000 A. Unrelated to this, strong 
Hq'-i-[NII] line emission is seen in both spectra and weak [SII] 
line emission in DC 303, in addition. Even though the observed 
excess in the TPN might be suspected to be caused by observa- 
tional effects or imperfect ISL SED modelling, it is not likely 
that at the same time DC 303 would show a much stronger ex- 
cess. If caused by problems in the ISL modelling, both glob- 
ules should show very similar excesses. We first summarise its 
observed properties and then discuss some alternative interpre- 
tations. 

The basic parameters characterising the emission feature 
are the central wavelength Ac-654?, A, width AA - 2800 A, 
band intensity /emission - 1-9 x 10"^ ergs cm"^ s"' sr"' for 
DC 303 and = 6950 A, = 1800 A, and /emission = 
4.7 X lO"*" ergs cm"^ s"' sr"' for the TPN. Compared to 
the scattered light within the emission band wavelength range, 
/emission //sea ~ 0.15 for DC 303 and 0.047 for DC 303 and the 
TPN, respectively. Emission band intensity and central wave- 
length were determined as the zeroth and first moments over 
the band; the width is an estimate made by eye. The emis- 
sion band intensity values are the observed ones. If corrected 
for Av = 0.5 mag foreground extinction, they would increase to 
2.7 and 0.67x10-^ ergs cm ^ s ' sr ' for DC 303 and the TPN, 
respectively. These values are similar to the estimate for the ERE 
band intensity in L 178 0, /ere ~ 2 x 10'^ e rgs cm"^ s"^ sr"' , 
as shown in Fig. 1 of Smith & WM (12002 ^ and, if interpreted 
in the same way as these authors, will lead to a similarly high 
photon conversion efficiency of ?; a: 10%. Our emission band 
to scattered light ratio for DC 303 is lower by factor of ~ 4 



than the estimate of ISmith & WiR d2002h for L 1780. Half of 
this difference is due to the stronger scattered light caused by 
the ~ 2 times stronger illumination from ISL in the direction 
1^1 = 15deg of DC 303 and the TPN as compared to ISL in the 
direction \b\ = 37deg of L 1780. 

The emission band intensities are, similarly as for the scat- 
tered light continuum, differences between globule rim and dif- 
fuse emission from the surrounding sky. Our observations do 
not allow us to put Umits to the emission band intensity at the 
sky positions. In their analysis of the scattered light spectrum of 
diffuse interstellar dust, Brandt & Draine (2012) have found no 
evidence for ERE emission. 

Strong ERE has been detected with good statistical and sys- 
tematic accuracy in high surface brightness objects, such as re- 
flection nebulae. Among the more difficult ERE searches in faint 
surface brightness objects the present result has some advan- 
tages. It is based on subtracting a scattered light spectrum, which 
is more realistic than in the previous works. These were based 
either on a general templa te galaxy spectr um as proxy for the lo- 
cal ISRF spectrum ( Szo moru et alJ[l998h . o r on a few broad- or 
intermediate-ban d photometric points only jGordon et al.ll 19981 : 
IWitt et alJl2008h . The present seai'ch for ERE in the two glob- 
ules can complement these previous determinations because of 
its wide spectral coverage and sufficient resolution. Especially, 
we suggest that our result is observationally p referable to th e 
low -resolution spectrophotome try of L 1780 by Mattila' ("1979^), 
which has been interpreted bv lChlewicki & Lau reiis (1987) to 
show ERE emission and has later been used as the only repre- 
sentati ve of an ERE-bearing dark/translucent cloud so far ( Smith 
& Witt l2002l : lGodard & Dartoisll2010l) . 

It is phenomenologically suggestive to interpret the excess 
emission at 5200-8000 A as extended red emission rising from 
the photoluminescence process. However, the interstellar en- 
vironment and the specific physical circumstances under which 
these globules were observed appear to speak against this inter- 
pretation: 

(i) Observations of reflection nebulae with different types of il- 
luminating stars have indicated that far-UV photons at wave- 
lengths shortward of 1 180 A are needed for ERE exitation (Dar- 
bon, Perrin, & Sivan 1999: Witt & Schild, 19 85: Witt et al. 2006J: 

(ii) these photons can penetrate only a thin surface layer of the 
globules: 

(iii) therefore, a uniform ERE band surface brightness profile 
over the globule's face is predicted rather than a bright rim - 
dark core structu re as observed for DC 303 and the TPN ( Witt. 
Oliveri, & Schildll990!): 

(iv) because our sky positions with Ay ~ 1 mag are also optically 
thick at far-UV wavelengths, practically equal amounts of ERE 
emission are expected at sky and globule sight-lines, causing the 
differential measurement to deliver a null result independently 
of the actual strength of the ERE emission. 

However, it is possible to construct somewhat ad hoc a model 
for the dust distribution and UV illumination in the globule sur- 
roundings that would make the ERE interpretation of the excess 
emission in DC 303 look more plausible: 

(i) Assuming that DC 303 is located behind most of the general 
Ay ~ I mag extinction layer as observed in the sky position. In 
this situation our differential measurement technique would not 
nullify the ERE signal because an equally high foreground ERE 
intensity is present both at the sky and the globule line of sight. 

(ii) DC 303 would then be exposed freely to the far-UV radia- 
tion field escaping from the OB stars in the Galactic disk. This 
radiation is impinging on DC 303 from the north (see Fig. 1) and 
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Table 3. Scattered and in situ line emission derived from the TPN and DC 303 spectra. In columns 8-10 the numbers in parentheses are the 
corresponding scattered continuum intensities in cgs units. The [Nii] line strengths include / i6583 and /i6548 l ines a nd [Sii] includes A6716 and 
A67 31 lines. In the lower part of the table the line equivalent widths are given as observed bvlBrandt & Drain3 i2012t) for diffuse dust and by Witt 
et al. ( 120101) for the dark nebula L 1780. 
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Table 4. Scattered and in situ emission line ratios. The [ Nii]/Hj3 ratios for the in situ case were calculated by assuming a ratio Hq'/H/3=2.87 
corresponding to case B for T = 10000 K an d no extinction lOsterbrock & Ferlan dll2006|). The TNii] line strengths include /16583 and A654S lines 
and [S ii] includes A6716 and /16731 lines. In lMadsen. Revn olds. & Haffner ( 2009) and lWitt et"aD ( l2010h the values are for single lines /16583 and 
/i6716 and were modified here corresponding to the line ratios. 
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would be similar in streng th to the ISRP that e xcites the ERE de- 
tected in the diffuse ISM (iGordon et al.lll998h . 
(iii) The missing ERE emission towards the dark core of DC 303 
can be understood since the far-UV radiation field is asymmet- 
ric, coming mostly from the side (north) and to lesser extent from 
the observer's direction. In this case the dark core would receive 
less far-UV illumination than the rim and less than assumed in 
the globule model of Witt, Oliveri, & Schild ( 1990). 

Two observational effects expected in this ad hoc model can 
be noticed: (i) the Galactic plane (north) side of DC 303, ex- 
posed to stronger far-UV radiation, is seen to be brighter than the 
south side in blue and red images (cf Fig. [U; (ii) if the TPN and 
DC 303 were both embedded in the suiTounding Ay ~ I mag 
extinction layer the SEDs of their bright rim diffuse emission 
should be the same for both and equal to the SED of ISL. On 
the other hand, if DC 303 is behind the Ay ~ I mag extinc- 
tion layer and the TPN in front or inside it, DC 303 should 
be more strongly reddened than the TPN. This is in fact the 



case: the colour of the bright rim of DC 303 is redder than 
that for the TPN, corresponding to an effective extinction of 
Ay ~ 1 mag (see Figs. Q and [8] upper panels). We also es- 
timated the extinction in the bright rim areas of the TPN and 
DC 303 using 2MASS JHK photometry of background stars in 
the two globule areas. The areas were defined as having a red 
surface brightness on DSS red plates (see Fig. [U of > 0.7x the 
maximum. The resulting R-band extinctions, determined from 
~24 stars in each globule, were 2.4 ± 0.3 and 3.0 + 0.3 for the 
TPN and DC 303, respectively. These values are compatible 
with the model where the bright rim corresponds to an optical 
depth through the globule of Ar ~ 2 mag, and the foreground 
extinctions are Ar ~ 0.4 mag for the TPN (embedded in the 
Ar ~ 0.8 mag layer) and Ar ~ 1 mag for DC 303 (behind the 
Ar ~ 0.8 mag layer). 

Another alternative to explain the excess emission at 5200- 
8000 A is in terms of an increased albedo over just this wave- 
length range and more strongly in DC 303 than in the TPN. Al- 
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Fig. 9. Excess emission from the bright rim of the TPN (lower spec- 
trum) and DC 303 (upper spectrum). The spectrum for DC 303 was 
shifted upwards by 10 cgs. The Hq'+[Nii] line emission in both spectra 
and [S ii] line emission in DC 303 are unrelated to the ERE. 

though such an albedo behaviour has not been observed before, 
it is also true that there have been very few albedo determina- 
tions in dense clouds covering this wavelength range with suf- 
ficient spectral resolution. On the other hand, the observations 
of reflection nebulae with late A to M-type illuminating stars 
shou ld have revealed, if presen t, such an albedo enhancement 
(Dar bon. Perrin. & Sivaniri999h . For B and early A-type illu- 
minating stars the ERE emission in reflection nebulae dominates 
this wavelength region and would hide such an albedo enhance- 
ment. 



7. Summary and conclusions 

We have obtained long-slit spectra 3500-9000 A of the faint sur- 
face brightness of two nearby {d ~ 200 pc) globules associated 
with the Chamaeleon dark cloud complex. The much stronger 
airglow emission spectrum and the atmospheric absorption lines 
were succesfully climated over most of the wavelength range. 
Spectra were presented separately for the bright rim and dark 
core regions of the globules. Although of relatively low resolu- 
tion {R ~ 100 - 200), the spectra allowed us to derive a number 
of results: 

The spectra are mainly caused by scattered light from 
the ISRF in the local Milky Way, consisting of a continu- 
ous/absorption line component of the integrated light of stars. 
The SED of the TPN bright rim, corresponding to an optical 
depth of T a; 2, is almost identical to the synthetic ISRF spec- 
trum. Also, in all globule spectra the stronger Fraunhofer lines 
of the integrated starlight spectrum were detected, with the lim- 
itation imposed by the relatively low spectral resolution of our 
observations. 

In addition, emission lines of Ha-i-[Nn], H/3, and [Sii] were 
detected and were interpreted in terms of scattered light from 
bright HII regions dominating the emission line component of 
the all-sky ISRF, plus an in situ WIM component behind the 
globules. Combining the line spectra of the bright rim and dark 
core regions allowed us to separate the scattered and in situ line 
components. 

After subtracting a model spectrum from the scattered light, 
there remains an excess emission between ~ 5500 - 8000 A in 
DC 303 and more weakly in the TPN, which is reminiscent of the 



photoluminescence band detected in reflection nebulae and other 
objects and is called ERE. The nature of this excess emission in 
DC 303 remains inconclusive however. 
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Appendix A: Synthetic spectrum of the integrated 
starlight 

Spectral synthesis is a common method in investigations of stel- 
lar populations in external galaxies that are too distant to be re- 
solved into individual stars. In this section we address the oppo- 
site problem: given the number densities and spatial distribution 
of the different types of stars in the solar neighbourhood, what is 
the spectrum of the integrated starlight (ISL) in different direc- 
tions of sky and for different vantage points of an observer (or 
a globule) off the Galactic plane ? In addition to the ISL, the 
Galactic surface brightness also contains the diffusely scattered 
starlight component, called the diffuse Galactic light (DGL). It 
contributes ca. 10-30% of the total light and has to be explicitly 
considered when estimating the total intensity of Galactic light. 
Because its spectrum is a copy of the ISL spectrum, it does not 
influence the strengths of the spectral features (absorption lines, 
bands, discontinuities), which are the main topic of our compar- 
ison of the ISL with the observed globule spectra. Our calcu- 
lation of the ISL spectru m between 300 and 10000 A follows 
the methods presented in lMattilal ( Il980allbh . We here give a brief 
description of ISL spectral modelling appropriate for the com- 
parison with the globule spectra. A more elaborate description 
and discussion of the results will be published in a separate pa- 
per. 



Appendix A.1: Galaxy model 

A simple model of the Galactic structure is adopted in which the 
stars and the dust are distributed in plane-parallel layers. The 
stars are divided into representative spectral groups and the num- 
ber density of stars per pc^, Z),(z), and the stellar emission coeffi- 
cient per pc^*, jiiz), in each group / are described by exponential 
density distributions in the z-direction: 



Diiz) = Di(Qy'iechiz/h,j) 



jiiz) = jiiO)sech(z/h^,i), 



(A.l) 



(A.2) 



where Z),(0) and ji(0) are the densities in the galactic plane 
and h^ i the scale height. The z-dependence of the extinction 
coefficient a{z) (mag kpc"') is approximated with the formula 



aiz) = a{Q)exp(\z\/fid), 



(A3) 



where we adopted a(Q) - 1.15 mag/kpc (in V band) and 
y8rf=134 pc. The extinction coefficient in the Galactic plane a{Q) 
is the average value for the distance range d = 0-2.5 kpc derived 
by Scheffler ( 1967a b) from a statistical analysis of the colour ex- 
cesses of 4700 stars near the Galactic plane (|z| < 75 pc) as pro- 
vided bylp eckel ( 1966). It a grees with the later estimates of e.g. 
iHakkil a.Mvers. & Stidham (1997) (a: 1 .5 maa/kpc) and Mendez 
& van Altena ( 1998). The scale hei ght was adopted from Mar- 
shall et al. (2006) and is in accordance (within ~ +20 pc) with 
other determinations for the solar neighbourhood within a few 
kpc. 

The results of lNeckell d 19661) . which we used in our 1980 
ISL modelling, have been derived for the more nearby Galactic 
environment, d < 500 pc, and indicated an extinction coefficient 
fl(0) = 2.46 mag kpc"' (V) and a scale height corresponding 
to /3d=4-0 pc. The ISL is, however, estimated to originate in a 
substantially higher volume around the Sun, reaching up to li > 
3 kpc in the Galactic plane. 

For the wavel ength dependence of intersteUar exti nction the 
formula given bv lCardelli. Clayton. & MadiisI ([19891) for Ry = 
3.1 was adopted. 

The dumpiness of the dust distribution has a profound effect 
on the ISL intensity. For a given mean extinction coefficient 
the mean ISL intensity for the clumpy dust distribution is higher 
than that for the homogeneous model. For a model consisting 
of a unifo rm distribution of stars and dust clouds Chandrasekhar 
& Miinch ([1950^ derived the following expression for the mean 
ISL intensity: 



1 



V 1 



-[1 



,-n(I-9) 



(A.4) 



Here v is the average number of clouds per kpc, n the mean 
number of clouds along the whole line of sight and q the average 
transmissivity of a dust cloud q = {exp{-T)), where t is the op- 
tical depth of one cloud. The parameters for the dust dumpiness 
correction were adopted according to the two-component model 
of lScheffleii (|1967a.b): there are two kinds of clouds with a mean 
optical thickness of f\{V) - 0.24 and f2{V) -1.5 and a relative 
fequency of vi - 0.9 and V2 = 0. 1 . 

The numerical integration along each line of sight was per- 
formed in steps of 10 pc and extended to Zmax-^6QQ pc, which is 
well outside the dust distribution layer in the directions \b\ + 0. 
For distances z > 3600 pc the ISL contribution was estimated 
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and found to be negligible. For b - deg the maximun distance 
along the plane was set to 10 kpc because high foreground ex- 
tinction (17.5 mag at z = 0, 8.3 mag at z = 100) prevents any 
significant starlight to penetrate from beyond this distance. For 
directions close to the Galactic plane, b < 10 deg, the maxi- 
mum distance to which the numerical integration was extended 
from the plane was limited to z„„„ <3600 pc by the condition 
not to exceed the maximun distance of 10 kpc along th e plane. 
For det ails of the model and the other formulae used see lMattilal 
(Il980ah . 



because the extinction coefficient at z = 50 pc is only 69% of 
its value at z = 0. The mean ISL intensity is seen to decrease 
modestly between z = and z = 50 pc. 



Appendix A.2: Stellar distribution parameters 

The division of stars into representative spectral groups cover- 
ing the different pa rts of the HR -diagram was made according 
to the approach of iFlvnn et al.l (12006) based on the analysis 
of the Hipparcos data base. Chris Flynn and Laura Portinari 
kindly provided us with their results in detailed tabular form. In 
these tables the stars are divided into the the following seven 
categories: (i, ii) main sequence (thin and thick disk), iv) 
clump stars (thin and thick disk), ( v, vi) old giants (thin and thick 
disk), and (vii) young giants. For each category the tables give 
as a function of My the stellar number density in the Galactic 
plane, D(0), and the scale height for a distribution of the form 
D{z) - D(Q)sech{z/h^). In addition, the colour indices B - V 
and V - Ic based on the Hipparcos and Tycho catalogues were 
derived for each magnitude slot of the seven categories. Because 
of the limited distance range of Hipparcos, its coverage for the 
supergiants was spa rse. We compleme nted this part by using the 
compilation of Wai nscoat et alJ (119921) . Our adopted parameters 
for the stellar distributions are shown in Table lA.4l 



Appendix A.3: Stellar spectral library 

For the comparison with the 3500-9000 A scattered light spectra 
of the globules we need a spectral library with good spectral type 
and wide wavelength coverage but with on l y a rel atively modest 
spectral resolution. The classical iPicklesI (Il998h library fulfils 
these conditions best. In order to choose the best template stars 
from the library for each of our stellar slots as given in Table lA.4l 
we used both the absolute magnitudes My and the colour indices 
B-V and V - as selection criteria. In most cases a satisfactory 
match with a Pickles template spectrum was possible. In a few 
cases we adopted the average of two nearby spectral templates. 
The Pickles catalogue numbers and spectral types for the chosen 
stellar templates are given in Table IA.41 



Appendix A.4: Results 

We show in Fig. lA. li the resulting ISL spectra for several galac- 
tic latitudes from \b\-0 to 90 deg for an observer located in the 
plane (z = pc). In Fig. lA.ll we also show ISL spectra for 
z = 50 pc, corresponding to the estimated distance of the two 
globules. The two spectra at b-ll.5 and 23.8 deg are bracketing 
the latitude range of the two globules, b=-14.2 deg for DC 303 
and b=-15.9 deg for the TPN. We also show the mean spectrum 
over the whole sky for the two z values. It can be seen that the 
absorption lines, bands, and discontinuities are very similar in 
all these ISL spectra. The ISL intensity for intermediate and 
high positive Galactic latitudes decreases, as expected, when the 
observer moves to z = 50 pc. On the opposite side of the sky 
(b < 0), not shown in the figure, the ISL intensity increases cor- 
respondingly. The ISL intensity at Galactic equator increases. 
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Dispersion axes 




Fig. 4. Example of a single spectrum of TPN as measured in the red arm of EMMI. Panel a shows the spectrum of diffuse light after background 
subtraction. The bright rim-dark core structure of the globule is evident. Panel b shows the total measured signal, i.e. the measured spectrum 
before the background subtraction when the dominating component is atmospheric emission. Note the positive correlation between the noise of 
the scattered light and the strength of the atmospheric emission. The extents of the bright rims and the dark core, as defined in Sect. 3.2, are 
marked with white lines. The location of the deep telluric absorption line at ~7650 A is marked with a ®-sign. The horizontal region located at 
~4.3' on the spatial axes is the gap between the two CCD arrays. The gap is interpolated along columns (spatial axes). The intensity of the 5577 A 
[01] atmospheric emission line is off-scale. Some stars outside the globule are seen as yellow strikes. The sky north of TPN (upper side of panel 
a) is dimmer than the sky south of TPN (lower side in panel a) because of the cometary shape of the cloud. 
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Fig. 5. Spectrum of diffuse light for TPN. Intensity A/ is given relative to the background sky adjacent to the globule. Panel a shows the 
spectra for the dark core and the mean of the bright rims. Panel b shows the ratio of the rim and the core spectrum after fitting and subtracting the 
continuum of the spectra and scaling the spectra to unity. The wavelengths of spectral lines and bands are marked with dashed lines. The spectral 
resolution of the measurements (PSF) is shown. The atmospheric transmission between ~ 6000 - 9000 A and the total measured signal (including 
the atmospheric emission lines) are shown in panel a for reference, in arbitrary scale. 



Article number, page 16 ofl20l 



K. Lehdnen and K. Mattila : Spectroscopy of diffuse light in dust clouds 




4000 5000 6000 7000 8000 9000 
?i[A] 



Fig. 6. Spectrum of scattered light for DC 303.8-14.2. The details are the same as for Fig.[5] 
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Fig. A.l. Integrated starlight spectra for different Galactic latitudes and for two positions of the observer, z = and z = 50 pc. The upper 
panels show the wavelength range 3500-10000 A, the lower panels show the wavelength range 3500-5500 A. For z = pc the spectra are 
shown from bottom to top for b = 90,36.8, 11.5, and deg as a black line. For z = 50 pc the spectra are shown from bottom to top for 
b = 90, 36.8, 23.6, 1 1 .5, 5.7, and deg as a black line. The mean spectra over the sky are shown as a red line. All spectra have the same scale and 
zero point. 
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Table A.l. Data for stellar distribution parameters 



No. 


Mv 


My range 




K 


./■() 


B-V 


V -I 


Pickles 


Pickles 




















Sp Class 


No. 




mag 


ma 




stars/pc"* 


pc 


My - stars/pc-* 


mag 


mag 






Main 


sequence, thin disk 
















1 


-4.0 


-4.20 


-3.25 


0.46E-06 


56. 


0.13E-04 


-0.28 


-0.28 


BO/Bl V 


3,4 


2 


-2.45 


-3.20 


-1.85 


0.34E-05 


56. 


0.32E-04 


-0.23 


-0.23 


Bl V 


4 


3 


-1.2 


-1.80 


-0.75 


0.13E-04 


56. 


0.38E-04 


-0.13 


-0.16 


B5-7 V 


6 


4 


-0.25 


-0.70 


0.20 


0.44E-04 


57. 


0.53E-04 


-0.10 


-0.10 


B8 V 


7 


5 


0.65 


0.25 


0.95 


0.92E-04 


68. 


0.52E-04 


0.01 


0.01 


AO V 


9 


6 


1.3 


1.00 


1.60 


0.15E-03 


82. 


0.44E-04 


0.14 


0.14 


A5 V 


12 


7 


1.8 


1.65 


1.95 


O.llE-03 


92. 


0.21E-04 


0.26 


0.26 


A7 V 


13 


8 


2.1 


2.00 


2.30 


0.14E-03 


107. 


0.20E-04 


0.32 


0.32 


A7/F0 V 


13,14 


9 


2.7 


2.35 


2.90 


0.41E-03 


140. 


0.36E-04 


0.42 


0.42 


FO V 


14 


10 


3.1 


2.95 


3.30 


0.42E-03 


176. 


0.24E-04 


0.49 


0.49 


F5 V 


16 


11 


3.5 


3.35 


3.75 


0.59E-03 


206. 


0.23E-04 


0.55 


0.56 


F6 V 


18 


12 


4.0 


3.80 


4.20 


0.71E-03 


220. 


0.18E-04 


0.59 


0.62 


F8 V 


20 


13 


4.4 


4.25 


4.55 


0.65E-03 


220. 


O.llE-04 


0.61 


0.66 


GOV 


23 


14 


4.7 


4.60 


4.90 


0.74E-03 


220. 


0.93E-05 


0.64 


0.69 


G2 V 


26 


15 


5.1 


4.95 


5.30 


0.92E-03 


220. 


0.82E-05 


0.70 


0.73 


G5 V 


27 


16 


5.5 


5.35 


5.70 


0.94E-03 


220. 


'S8F-0'S 


0.78 


0.81 


G8 V 


30 


17 


5.9 


5.75 


6.15 


O.llE-02 


220. 


0.45E-05 


0.88 


0.91 


KO V 


31 


18 


6.4 


6.20 


6.90 


0.19E-02 


220. 


0.47E-05 


0.94 


0.99 


K2 V 


33 


19 


7.35 


6.95 


8.05 


0.33E-02 


220. 


0.35E-05 


1.10 


1.27 


K4 V 


35 


20 


8.8 


8.10 


9.35 


0.46E-02 


220. 


0.15E-05 


1.42 


1.77 


MOV 


38 


21 


9.9 


9.40 


11.10 


O.llE-01 


220. 


0.94E-06 


1.50 


2.17 


M2.5 V 


41 


22 


12.2 


11.15 


12.95 


0.15E-01 


220. 


0.27E-06 


1.62 


2.64 


M4 V 


43 


Main sequence, thick disk 
















23 


2.1 


2.00 


2.30 


0.39E-06 


670. 


0.50E-07 


0.27 


0.27 


A7 V 


13 


24 


2.7 


2.35 


2.90 


0.45E-05 


670. 


0.38E-06 


0.37 


0.37 


FO V 


14 


25 


3.1 


2.95 


3.30 


0.80E-05 


670. 


0.44E-06 


0.44 


0.44 


F2 V 


15 


26 


3.5 


3.35 


3.75 


0.50E-04 


670. 


0.18E-05 


0.50 


0.51 


wF5 V 


17 


27 


4.0 


3.80 


4.20 


0.78E-04 


670. 


0.20E-05 


0.54 


0.57 


wF8 V 


21 


28 


4.4 


4.25 


4.55 


0.72E-04 


670. 


0.13E-05 


0.56 


0.61 


F8 V 


20 


29 


4.7 


4.60 


4.90 


0.82E-04 


670. 


O.lOE-05 


0.59 


0.63 


GOV 


23 


30 


5.1 


4.95 


5.30 


O.lOE-03 


670. 


0.91E-06 


0.65 


0.68 


wG5 V 


28 


31 


5.5 


5.35 


5.70 


lOF-03 

V/. J. \Jl—i \JJ 


670. 


6'SF-06 

\J .\J^ l—i \J\J 


0.73 


0.76 


G8 V 


30 


32 


5.9 


5.75 


6.15 


0.12E-03 


670. 


0.50E-06 


0.83 


0.86 


KO V 


31 


33 


6.4 


6.20 


6.90 


0.21E-03 


670. 


0.52E-06 


0.89 


0.94 


K2 V 


33 


34 


7.35 


6.95 


8.05 


0.37E-03 


670. 


0.38E-06 


1.06 


1.22 


K4 V 


35 


35 


8.8 


8.10 


9.35 


0.51E-03 


670. 


0.17E-06 


1.37 


1.72 


MOV 


38 


36 


9.9 


9.40 


11.10 


0.13E-02 


670. 


O.lOE-06 


1.45 


2.12 


M2/2.5 V 


40,41 


37 


12.2 


11.15 


12.95 


0.17E-02 


670. 


0.30E-07 


1.56 


2.59 


M3/M4 V 


42,43 


Clump stars, 


thin disk 


















38 


-1.6 


-2.10 


-1.15 


0.36E-06 


180. 


0.13E-05 


2.1 


2.1 


M3 III 


98,113 


39 


-1.0 


-1.10 


-0.95 


0.25E-06 


180. 


0.63E-06 


1.77 


1.76 


Ml III 


96 


40 


-0.85 


-0.90 


-0.80 


22F-06 


180. 


49F-06 


1.57 


1.60 


KS/MO TIT 


93,95 


41 


-0.7 


-0.75 


-0.60 


0.34E-06 


180. 


0.64E-06 


1.45 


1.45 


K5III 


93 


42 


-0.5 


-0.55 


-0.45 


0.27E-06 


180. 


0.43E-06 


1.26 


1.26 


wK4 III 


91 


43 


-0.35 


-0.40 


-0.30 


0.27E-06 


180. 


0.37E-06 


1.16 


1.14 


wK3 III 


88 


44 


-0.2 


-0.25 


-0.10 


0.52E-06 


180. 


0.61E-06 


1.02 


1.02 


K2III 


85 


45 


0.3 


-0.05 


0.65 


0.29E-04 


180. 


0.18E-04 


0.88 


0.89 


G5II/WG8III 


110,77 


46 


0.95 


0.70 


1.20 


0.52E-04 


180. 


0.25E-04 


1.05 


1.05 


wK2 III 


85 


Clump stars. 


thick disk 


















47 


-0.7 


-0.75 


-0.60 


0.44E-08 


670. 


0.78E-08 


1.45 


1.45 


K5III 


93 


48 


-0.5 


-0.55 


-0.45 


0.60E-07 


670. 


0.94E-07 


1.26 


1.26 


wK4 III 


91 


49 


-0.35 


-0.40 


-0.30 


0.18E-06 


670. 


0.25E-06 


1.15 


1.14 


wK3 III 


88 


50 


0.4 


-0.25 


1.00 


0.12E-04 


670. 


0.65E-05 


0.90 


0.95 


WG8III/G5II 


77,110 
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No. 


My My range 


Do 




;■() 


B - V 


V - 1 


Pickles 


Pickles 


















Sp Class 


No. 




mag mag 




stars/pc 


pc 


My - Stars/pc-' 


mag 


mag 






Oldg 


;iants, thin disk 


















51 


-0.6 -1.25 


0.00 


0.18E-04 


180. 


0.29E-04 


1.45 


1.45 


K5III 


93 


52 


0.5 0.05 


1.00 


0.24E-04 


180. 


0.16E-04 


1.31 


1.28 


wK4 III 


91 


53 


1.5 1.05 


2.00 


0.50E-04 


180. 


0.12E-04 


1.07 


1.06 


Kl III 


81 


54 


2.5 2.05 


3.00 


0.67E-04 


180. 


0.65E-05 


0.95 


0.95 


KOIV 


57 


55 


3.5 3.05 


4.00 


0.13E-03 


180. 


0.49E-05 


0.86 


0.86 


G8 IV 


56 


56 


4 15 4 05 


4.65 


0.64E-05 


180. 


0.15E-06 


0.85 


0.85 


G8 IV 


56 


Oldg 


;iants, thick disk 


















57 


-1.0 -2.00 


0.00 


0.21E-05 


670. 


0.50E-05 


1.45 


1.45 


K5 III 


93 


58 


0.5 0.05 


1.00 


0.24E-05 


670. 


0.14E-05 


1.28 


1.28 


K4 III 


90 


59 


1.5 1.05 


2.00 


0.65E-05 


670. 


0.15E-05 


1.07 


1.06 


wKl III 


82 


60 


2.5 2.05 


3.00 


0.15E-04 


670. 


0.15E-05 


0.95 


0.95 


KO IV 


57 


61 


3.25 3.05 


3.45 


0.99E-05 


670. 


0.49E-06 


0.88 


0.88 


G8 IV 


56 


Youn 


g giants, thin disk^ 




















-0.5 -1.00 


0.00 










1 44 


KA TTT 

XVT- 111 




63 


0.5 0.05 


0.95 


0.57E-05 


56. 


0.38E-05 


0.87 


1.28 


K3 III 


87 


64 


1.7 1.00 


2.45 


0.84E-05 


100. 


0.17E-05 




1.03 


Kl III 


81 


65^ 


-5.14 




0.20E-06 


56. 


2.28E-05 






05/B5 V 


1,117 


66"= 


-6.05 




0.33E-07 


56. 


0.84E-05 






A2/G0/G2 I 


120, 122, 125 


67'' 


-6.34 




0.32E-07 


56. 


l.lOE-05 






K2/M2I 


128,131 


68*= 


-5.79 




0.13E-07 


56. 


0.27E-05 






M3 II 


113 



Notes. Parameters for supergiants are according to Wainscoat et al. (1992) OB A-G I-II K-M2 I-II M3-M4 I-II 
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